v±: A small, but sensitive observable for heavy ion collisions 
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Higher order Fourier coefficients of the azimuthally dependent single particle spectra resulting from 
noncentral heavy ion collisions are investigated. For intermediate to large transverse momenta, these 
anisotropics are expected to become as large as 5 %, and should be clearly measurable. The physics 
content of these observables is discussed from two different extreme but complementary viewpoints, 
hydrodynamics and the geometric limit with extreme energy loss. 
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Collisions of large nuclei at nonzero impact parameter 
exhibit the special feature of a strongly deformed over- 
lap region. The subsequent dynamical evolution of the 
collision zone converts the initial spatial eccentricity into 
an azimuthal anisotropy of the final state observables. 
In particular anisotropies in momentum space [1] have 
been analyzed experimentally in great detail in recent 
years, as they have been shown to be generated dur- 
ing the earliest and hottest stages of the reaction [2]. 
As momentum anisotropies are generated, the eccentric- 
ity in coordinate space is smoothed out, bringing fur- 
ther generation of momentum anisotropies to a stall, a 
process that occurs within about the first 5 fm/c of the 
total 15 fm/c lifetime of the fireball [2-4]. In Au+Au 
collisions at RHIC (the Relativistic Heavy Ion Collider 
at Brookhaven National Laboratory), those anisotropies 
have been found [5] to be as large as predicted by hydro- 
dynamic calculations that assume a short equilibration 
time T C q U < 1 fm/c [4,6]. Further theoretical investiga- 
tion has revealed the sensitivity of anisotropic flow on the 
equation of state of the expanding medium [7-9] , indicat- 
ing the necessity of a hard phase at high temperatures 
with a soft transition region of width Ae ~ 1 GeV/fm 3 to 
lower temperatures in order to describe the data in more 
detail. For these reasons anisotropic flow became one of 
the most attractive tools to study the nuclear equation of 
state of the matter created in the collision by utilizing the 
most abundant particles, alternatively to study statisti- 
cally disadvantaged rare probes with which one tries to 
infer the properties of the bulk by modeling their mutual 
interaction with the medium they traverse. Other predic- 
tions from hydrodynamic concepts such as centrality and 
mass dependence of momentum anisotropies [8] have sub- 
sequently been confirmed to be in qualitative and quan- 
titative agreement with experiments [10,11], as long as 
the applicability of hydrodynamics is not overstretched, 
i.e. for impact parameters not larger than b ~ 7 f m and 
transverse momenta not to exceed px ~ 1.5 Ge V. Beyond 
these values clear deviations from the hydrodynamic pre- 
dictions start to occur. This is expected as the smaller 
systems do not sufficiently equilibrate and the most rapid 
particles escape the fireball without fully participating in 
the collective dynamics. 

To date, anisotropies at midrapidity are generally only 
characterized in terms of the second coefficient v-i of the 
Fourier decomposition of the azimuthally sensitive mo- 



mentum spectrum. With the z direction given by the 
beam axis and the x direction defined by the direction 
of the impact parameter (and y perpendicular to both), 
the general expression for the Fourier decomposition with 
respect to the azimuthal angle ip — arctan(y/x) is 
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1 + ^2 2w„(pT)cos(n^) , (1) 
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where no sine terms appear due to the symmetry with re- 
spect to the reaction plane, and all odd harmonics vanish 
at midrapidity due to the symmetry ip <-> ip + ir. (Note 
that the latter is not true for noncentral c?-Au collisions, 
where v\, V3,... should be useful quantities to classify 
and study the measured spectra - even at midrapidity). 
In the following we will elaborate on the centrality and 
momentum dependence of higher anisotropy coefficients, 
which turn out to exhibit interesting features that, with 
the currently available data from RHIC's run in 2002, 
might already be accessible. 

Generally, higher momentum anisotropies are expected 
to be small, and in fact no experimental data have been 
published until this day. An early hydrodynamic cal- 
culation [12] predicted the momentum integrated value 
of V4 not to exceed 0.5 % at SPS energies even in the 
most peripheral collisions, and shows only a slight in- 
crease with beam energy up to RHIC energies. In another 
study Teaney and Shuryak [13] calculated a±, which is 
the fourth coefficient of the momentum integrated differ- 
ential particle spectra weighted by an additional momen- 
tum squared (which was applied to enhance the effects 
of transverse flow). This quantity is thus strongly biased 
to larger px- From their hydrodynamic calculation they 
obtain oti ~ 1%; clearly the momentum integrated value 
of i>4 is still smaller. 

Today, experiments at RHIC achieve sufficient statis- 
tics to measure the elliptic flow coefficient differentially 
in momentum out to large transverse momenta up to 10 
GeV [14]. In midcentral collisions this coefficient is found 
to saturate starting at px ~ 2 GeV at a large value of 
20 - 25%, indicating that three times as many particles 
are emitted into the reaction plane than out of the reac- 
tion plane. Such a large value also means that the first 
order deformation of Eq. (1) is not elliptic anymore, but 
that a polar plot of the azimuthal distribution resembles 
more a peanut than an elliptical shape, aligned with its 



1 



longer axis into the reaction plane, see Fig. 1. With- 
out any higher order coefficients, the elliptic deformation 
starts to develop such a waist when v-i becomes larger 
than 10 %, which can be shown explicitly by discussing 
the term l+cos2(p in Cartesian coordinates (x, y) around 
i = 0. 
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FIG. 1. (Color online) Polar plot of the azimuthal distri- 
bution of different Fourier expansions. Left, distortion of the 
unit circle through an elliptic component V2, shown for 10% 
(green) and 20% (red). Right, fourth component Va for the 
same values. 

This raises the interesting question, whether at high 
Pt Fourier coefficients of higher order might become suf- 
ficiently large to restore the elliptically deformed shape 
of the particle distribution. To restore the elliptic shape 
as much as possible only through introduction of a fourth 
anisotropy coefficient, V4 ~ gj(10 v 2 — 1) is required, as 
is derived straightforwardly from Eq. (1).* This, in con- 
sideration of the large experimental values for v%, might 
indicate that 1)4 could reach values around 3-5 %. 

To investigate higher flow coefficients in more detail, 
we first revisit hydrodynamic calculations, although they 
find their limitations for pions of transverse momenta 
beyond 2 GeV, and for heavier baryons at momenta 3- 
4 GeV (for a recent review see Ref. [15]). Later on we 
will discuss a complementary model in terms of a sim- 
ple picture for the limit of very large transverse mo- 
mentum. The results presented in the following were 
achieved by a more modern hydrodynamic calculation 
than in Ref. [12], as it becomes important at RHIC ener- 
gies to correctly account for features of chemical non- 
equilibrium during the hadronic evolution of the sys- 
tem [16,17], a feature that significantly influences the 
relation of energy-density and temperature [18,19]. At 
temperatures beyond T cr ;t = 165 MeV the underlying 
equation of state turns into a hard ideal gas equation of 
state via a strong first order phase transition in order to 
mimic the transition from a gas of interacting resonances 



*The definition of restored ellipticity is ambiguous. Here we 
chose to find the value of V4 where the second derivative of 
y(x) vanishes at x = 0, thus giving rise to a very smooth 
transition across ip — tt/2. 



to a hard plasma phase. The initial conditions of this 
calculation were determined to fit single particle spectra 
measured in Au+Au collisions at ^/snn = 200 GeV. For 
further details, see Ref. [17]. 

Figure 2 shows the momentum dependence of the 
Fourier coefficients of pion spectra up to order 8. "El- 
liptic flow," i>2 dominates the emission anisotropy at all 
momenta. (To use the same vertical scale for all coeffi- 
cients, the second harmonic V2 was divided by 10.) As an- 
ticipated, V4 increases with momentum, reaching a value 
of about 2.8 % at px — 2 GeV (thus not quite enough to 
restore the particle distribution to an elliptic shape, as 
V2(pt = 2 GeV) w 29%). Also still has a sizable nega- 
tive value of about 1.2%. The mass dependence of higher 
flow harmonics does not bring any surprises. Due to their 
larger masses, heavy particles pick up larger transverse 
momenta when they participate in the established collec- 
tive flow. As for the spectra and elliptic flow, also the 
higher anisotropies are shifted out toward larger trans- 
verse momenta. Thus at a given transverse momentum 
below 2 GeV, the anisotropies of heavy resonances are 
smaller than of pions. Comparing hydrodynamic calcu- 
lations to experimental results, it appears however that 
deviations from hydrodynamic behavior occur later for 
heavy particles, and therefore anisotropies of baryons can 
eventually surpass anisotropies of light mesons [20] . From 
these calculations one can expect a fourth flow compo- 
nent of (anti-)protons of about 7 % at 4 GeV transverse 
momentum, a sizable quantity that should be measur- 
able. 
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FIG. 2. (Color online) Fourier coefficients of the pion spec- 
tra resulting from a hydrodynamic calculation to describe 
particle spectra at ^/snn = 200 GeV at impact parameter 
6 = 7 fm. «2 is scaled by a factor 0.1. 

Surprisingly it turns out that 114 is highly sensitive to 
the initial conditions of the calculations: initializing the 
calculation with a prehydrodynamic, isotropic radial flow 
field as in Ref. [17] does not change the magnitude of the 
flow anisotropies by a large amount. However we find 
that it completely changes the sign of V4, while the other 
coefficients remain largely the same. Clearly W4 at low 
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and intermediate pr has therefore a strong potential to 
constrain model calculations and carries valuable infor- 
mation on the dynamical evolution of the system. To 
study how the higher flow anisotropies develop in the 
course of the hydrodynamic evolution we calculate mo- 
ments of the transverse flow field (v x , v y ) through 



S n (r) = 



J dxdy e(x, y; r) 7 vt cos(nip v ) 
J dxdye(x,y;r)-f 



(2) 



where e{x,y 1 r) is the energy density in the transverse 
plane at time r and = {v*. + v 2 ) 1 / 2 is the transverse 
flow velocity at a given point (x,y), 7 = (f — v\)^ 1 l 2 1 
and <p v — &rct&n(v y / v x ) is the angle of the local flow 
velocity with respect to the reaction plane. (Note that 
for n = this reduces to the definition of the mean ra- 
dial velocity used in earlier studies [4], but for n = 2 
it is slightly different from the momentum anisotropy e p 
which is defined in terms of the difference in the diag- 
onal elements of the energy momentum tensor.) Figure 
3 displays the time evolution of the flow coefficients and 
clearly shows that radial flow Sq increases throughout the 
lifetime of the system (which is about 15 fm/c), with a 
somewhat reduced acceleration while most of the system 
is in the mixed phase at around 4 fm/c. Quite oppositely 
anisotropy components quickly saturate and remain con- 
stant, as the initial geometric deformation of the source is 
lost. About 2/3 of the final value of 84 is generated dur- 
ing the earliest stages of the evolution, where the system 
is governed by the hard equation of state of the quark 
gluon plasma. Sq quickly achieves small negative values 
of the order of -0.2 % (as expected from the negative value 
of vq). Due to its small values, it is subject to numeri- 
cal fluctuations and therefore not shown. It is expected 
that the flow coefficients are monotonically related to the 
anisotropy coefficients in the final particle spectra. De- 
tails of this relation are subject to future study. 
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FIG. 3. (Color online) Time evolution of the anisotropy 
coefficients in the transverse velocity field as given by the 
hydrodynamic calculation. 5q characterizes the radial flow, 62 
the elliptic deformation and 84 the quadrupole deformation. 

In a model complementary to the hydrodynamic 
description which is constrained to the low momen- 



tum, thermalized part of the spectrum, one can study 
anisotropies that result from opacity effects that the de- 
formed initial geometry exerts on particles of high mo- 
menta [21-23]. Assuming extreme opacity, or extreme 
jet quenching, particles of large transverse momenta from 
any part of the surface can only be emitted within ±90° 
of the normal vector to the surface, which for this ex- 
ploratory study is assumed to be given by the overlap of 
two hard spheres (radius R)J In this picture, the nth 
Fourier coefficient of the resulting expansion is found to 
be (for n even) 



1 — n z na 



(3) 



where a = arccos(6/2i?) characterizes the centrality in 
terms of the opening angle between the line through the 
center of one nucleus to the intersection point of the two 
nuclear spheres and the reaction plane. Note that as 
b increases from to 2R, a decreases from ir/2 to 0. 
For n = 2, V2 — sin(2a)/(6o;), a formula first given by 
Voloshin [26]. For n > 2, Eq. 3 becomes significantly 
more interesting (unfortunately, however, the size of the 
expected signal drops with n~ 3 ). Note first the oscillat- 
ing sign of the prefactor which will manifest itself in an 
oscillation of the coefficients at large b (small a). Sec- 
ond, one notices that for larger n the sine function starts 
changing sign upon scanning impact parameters (b <-> a), 
whereas for n — 2, the Fourier coefficient increases lin- 
early with impact parameter from to 1/3, and is thus 
always positive. In particular for n — 4 one observes a 
transition from positive to negative values at b — V%R, 
reaching the limiting value V4 — — 1/15 « —6.67% for 
b — > 2R, again a rather large number. The local max- 
imum of V4 is at b 0.865i? with V4 « 1.45 %. The 
centrality dependences of the second to eighth Fourier 
coefficients are shown in Fig. 4 (as before vi is divided 
by a factor of 10). 

Clearly, the expected values for higher expansion coef- 
ficients from this calculation are small. Applying a more 
realistic model does, however, not necessarily lead to a 
further reduction of the signal, but might in contrast in- 
crease, similarly as observed for calculations of elliptic 
flow [24] . This is also suggested by the experimental data 
on V2 which seem to exceed the limits given even by the 
most optimistic assumptions for anisotropies from geom- 
etry and energy loss. In any case, such a distinctive signal 
as a sign change upon variation of centrality should be 
preserved when model calculations get refined and should 
be accessible within the statistics reached by today's ex- 
periments. 



1 Relaxing this geometrical constraint by adopting a finite 
diffuseness of the nuclei reduces the resulting asymmetries, far 
below the experimental results, and novel effects are required 
to explain those large limiting values [24,25]. 
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FIG. 4. (Color online) Centrality dependence of the 
anisotropy coefficients in the limit of pure surface emission 
from the deformed overlap region to describe the saturation 
value of the anisotropies at large transverse momenta. V2 is 
scaled by a factor of 0.1. 

As an aside we would like to discuss a new but naive 
experimental way of measuring higher order flow coef- 
ficients, in particular V4. Apart from the "straightfor- 
ward" method to determine each event's reaction plane 
and averaging all particles' cos(4<£>) with respect to the 
reaction plane and the more recently developed cumu- 
lant technique [27], simply measuring the particle spec- 
tra with respect to the reaction plane might be sufficient 
(theoretically). Having the azimuthally sensitive particle 
spectra given, Eq. (1) shows directly that [abbreviating 
dN/p T dp T dyd(p(p T , <p) = n(<p)] 

n(0) + n(7r/2) - 2n(7r/4) = 

= T — C ^—[&V4(pt) + 8v 1 2(pt) + --}- (4) 
27T pxdpx 

From the experience gathered so far we expect U12 to be 
negligibly small compared to U4. Finite opening angles 
for the azimuthally sensitive distribution will reduce the 
signal again, however, the V4 is in this approach enhanced 
by a factor of 8, which might still render it measurable. 
We propose this method thus as an interesting indepen- 
dent approach to measure V4 (but also , by simply sub- 
tracting spectra in and out of the reaction plane). 

In summary, we have shown that v^{pt) might achieve 
significantly large values at intermediate and large trans- 
verse momenta in noncentral heavy ion collisions and 
contains important physics of heavy ion collisions, a fact 
largely overlooked in the past. Numerical estimates for 
its size were given both in the limit of full thermaliza- 
tion with a subsequent hydrodynamic expansion as well 
as in the geometrical limit of extreme "jet quenching." 
For the latter case a strong centrality dependence and a 
change of the sign of the signal for higher anisotropies 
was found. From both approximations we found that V4 
might reach values of the order of 5%, clearly observ- 
able with todays equipment and statistics. 1*4 has been 
found to be highly dependent on the initial configuration 



of the system and its evolution and is thus an important 
new tool to constrain model calculations and analyze de- 
tails of the system's history. Consequently it deserves to 
be studied in great detail both experimentally as well as 
theoretically. 
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